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Abstract 



This study elaborates some examples of a simple evolutionary stochastic 
rate process where the population rate of change depends on the distribution 
of properties — so different cohorts change at different rates. We investigate the 
effect on the evolution arising from parametrized perturbations of uniformity 
for the initial inhomogeneity. The information geometric neighbourhood system 
yields also solutions for a wide range of other initial inhomogeneity distributions, 
including approximations to truncated Gaussians of arbitrarily small variance 
and distributions with pronounced extreme values. It is found that, under quite 
considerable alterations in the shape and variance of the initial distribution of 
inhomogeneity in unfitness, the decline of the mean does change markedly with 
the variation in starting conditions, but the net population evolution seems 
surprisingly stable. 
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1 Introduction 

Consider a population with an inhomogeneous property distribution faced with a 
sudden environmental change at t = and let N(t) represent the declining population 
of more unfit individuals. In our example the distribution of unfitness a will lie in the 
range [0, 1] and under selective evolution the fraction N(t) could consist of individuals 
with unfitness below some threshold value which itself may evolve. However, as 
N(t) evolves we expect the distribution of unfitness to become skewed increasigly 
towards smaller values, so improving population fitness. We make use of existing 
inhomogeneous stochastic rate process theory |5] , and information geometry, [TJ [21 
IH [5] . In the example we elaborate, the inhomogeneous population N is classified by 
a smooth family of probability density functions {Pt,t > 0} with random variable 
< a < 1, having mean E t {a) and variance a 2 (a) = E t (a 2 ) — (E t (a)) 2 . Here a 
represents an unfitness that controls the decline in the frequency of the a-cohort, so 
we could view fitness as the variable 1 — a. 

Let lt[a) represent the frequency at the a-cohort, then we have 
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Figure 1: The log-gamma family of probability densities Po(a,v,T) from (11) as a 
surface for the case of central mean Eo(a) = h. This symmetric surface coincides 
with the uniform density 1 at r = 1, and tends to the delta function as r — > oo 



Karev [5] obtained general solutions for these equations giving us 

/>oo 

N{t) = N(0)L (t) where L (t) = / P (a) e - at da (3) 



dN , s f°° d\ogL 

-E t {a)N where E t (a) = / aP t (a)da = 1— ^ (4) 



eft ' In dt 



dE t (a) 
dt 



o-i{a) = {E t {a)y-E t {a*) (5) 



P t (a) = e~ atI ^l and l t (a) = e~ at L (t) (6) 



dt 



P t (a){E t (a) - a). (7) 



Here Lo(^) is the Laplace transform of the initial probability density function Pn(a) 
(which is of course zero outside [0, 1]) and so conversely Po(a) is the inverse Laplace 
transform of the population (monotonic) decay solution . See Feller [H| for more 
discussion of the existence and uniqueness properties of the correspondence between 
probability densities and their Laplace transforms. We see from (J7]) that when the 
unfitness a exceeds its mean E t (a) then the population density of that cohort declines 
and conversely the densities of cohorts with a < E t (a) tend to grow. The Shannon 
entropy at time t is 

St = -Ek (logP t (a)) = -Ek (log P0i ^ at ) (8) 

which reduces to 

S t = So + log L Q (t) + E t (a)t. (9) 
By using = — a\ (a), the decay rate is then 

d dT = - (j2{t)t - {10) 

This shows how the variance controls the entropy change during quite general inho- 
mogeneous population processes, as we saw in [5]. 

Karev gave the particular solutions for the cases of initial densities that were 
Poisson, gamma or uniform. In [5j we studied the case of a bivariate gamma density 
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Figure 2: Log-gamma probability density functions Po{a) from for a G [0, 1], with 
central mean E n (a) — ~, and r from 0.2 to 2 in steps of 0.2. Note that the parameter 
t controls the shape of the graph and for r = 1 we have Po(a) = 1. 




Figure 3: Mean and variance for the log-gamma family. 



for an epidemic situation. Here we shall study a family of log-gamma densities that 
determine a neighbourhood of the uniform distribution ;!2j , so recovering the solution 
of Karev for the uniform distribution as a special case [y = r = 1). The log-gamma 
family contains also close approximations to Gaussians of arbitrarily small variance 
truncated to [0, 1] (for r > 1) as well as densities having outlier cohorts of extreme 
values that are enhanced (r < 1), Figure [I] 

2 Evolution of inhomogeneity from an initial log- 
gamma probability density 

We studied in [2] the smooth family of log-gamma distributions and their information 
geometry. This family has probability density function 

a v - x v r |log(i)r _1 

for random variable a € [0,1] and parameters v, r > 0. The mean and variance, 
Figure [3j are given by 

E » ( "» = {ihj (12) 
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Figure 4: An affine immersion in M 3 of the 2-manifold of log-gamma probability 
densities. The black curves in the immersion represent the distributions with v = 1 
and t = 1 and centred on their intersection is a spherical neighbourhood of the uniform 
distribution. The other two points are for the two cases (v = 0.1, r = 0.289) and 
(v = 2.75, r = 2.24) cf. Figures |6|7|a 



The locus in this family of those with central mean Eo(a) = | satisfies 

i/(2* - 1) = 1 (14) 

and some are shown in Figure [TlThe uniform density is the special case with r = v = 1 
and this can be seen in Figure[2j with sections through the surface of Figure [T] The 
log-gamma family yields a smooth Ricmannian 2-manifold with coordinates {v, r) € 
R + x R + and metric tensor given by 

(15) 

In fact, this manifold is an isometric diffeomorph of the 2-manifold of gamma densities, 
with random variable x = —log a via natural coordinates {y, r) (cf. [1]) for which 
t = 1 corresponds to the subfamily of exponential densities. Through this smooth 
diffcomorphism we can therefore represent the manifold of log-gamma densities as a 
natural affine immersion (cf. [2]) of K + x K+ in M 3 

{v,t) t-> {i/,r,logr(r) -rlogi/}. (16) 

This is illustrated in Figure [4] which shows also a spherical neighbourhood in R 3 
centred on the point at the uniform density, v — 1 — r; the other two points are for 
the two cases (y = 0.1, t = 0.289) and (y = 2.75, r = 2.24) which represent densities 
also having mean value | and used in the sequel as initial unfitness distributions, cf. 
Figures [6|7pl 

The log-gamma entropy ^ is given by 

Slg(v,t) - i/ t (i/+1)-''- 1 (A(i/,t) + B(i/,t)) (17) 

with A(z/,r) = r(i/+(i/ + l)log(z/ + l)-l)-(i/+l)(r-l)V'(T) 

, / x / /^ _r r(r)\ x dlogr(r) 

andB( i /,r) = {v + 1) log ( - - ^ M , i>(r) = ^ ' . 

Figure [5] shows the graph as a surface and as a contour plot with gradient vector 
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Figure 5: Shannon entropy function for the log-gamma family as a surface (left) and 
as a contour plot with entropy gradient flow (right). 



field indicated; the main curvature occurs near the origin. 

The Laplace transform integral ([3]) for the log-gamma density with general (v, r) 
seems intractable so we used a series development up to sixth order for the term 
e~ at , bearing in mind that a G [0, 1]. Thus, up to sixth order in t > we obtain a 
good approximation to Lo(t), accurate to about 0.1% up to t = 1. The corresponding 
expressions for Lo(t), Pt(d), E t (a), of (a) and hence N(t) are known but somewhat 
cumbersome to present here so we present some graphics for illustration. The decline 
of the mean E t (a) as the density P t (a) develops from three initial log-gamma densities 
Po(a) with central mean, Eq(o) = | is shown in Figure [6] The more rapid decline 
occurs for the initial density with r = 0.289 and the slower decline occurs for the case 
having r = 2.24 with the uniform case r = v = 1 in between, cf. also Figure [2] The 
corresponding decline of N(t)/N(0) = Lo(t) is shown in Figure [7] and the variance is 
shown in Figure [8] 

We do have an analytic solution for the special 1-parameter family of log-gamma 
densities with r = 1 (the 'log-exponential' densities), from ^ and 

e~ at a u ~ 1 f 

= r(,)-r(M) (18) 

r(i/ + i) -T(v + i,t) , , 

Et{a ^ l) = Vm-Am)) (19) 



m 

N(0) 



Lo^^ut-^T^-T^t)) (20) 



where T(is, t) is the incomplete gamma function. E t (a, 1), declines for all choices of 
parameter v in the initial density, which includes of course also the uniform density. 
Further analytic solutions arise for other initial log-gamma densities with positive 
integer r (the 'log- Pearson Type III' distributions) cf. Figure [2j as generalized hyper- 
geomentric functions; for example, 

r = 2 =► =L (t)= 2 F 2 (i/,i/;i/ + l,i/ + l;-i) (21) 

r = 3 => W-=Lo(t)= 3 F 3 {v,v,v;v + 1, i/ + 1, i/ + 1; -t). (22) 
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Figure 6: Initial log-gamma densities Po(a) shown in left panel with central mean 
Eq(o) = \ for the cases r = 0.289, v = 0.1 (lower graph), r = 2.24, v = 2.75 (upper 
graph), and the uniform density Pq(o) = 1 for r = v = 1. TTie panel shows the 
decline with time of the mean E t (a) from for these initial densities. 
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Figure 7: Initial log-gamma densities Po(a) shown in left panel with central mean 
Eq(o) = \ for the cases r = 0.289, v = 0.1 (lower graph), t — 2.24, v = 2.75 (upper 
graph), and the uniform density Po(a) — 1 for r = v = 1. TTie ri^/ii panel shows the 
fractional decline with time of the population N(t)/N(0) = Lo(t) from £3|) for these 
initial densities. 
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Figure 8: Initial log-gamma densities Po(a) shown in left panel with central mean 
Eo(a) — | for the cases r = 0.289, v = 0.1 (lower graph), r = 2.24, v = 2.75 (upper 
graph), and the uniform density Pq(o) = 1 for t — v = 1. The right panel shows the 
decline with time of the variance of (a) from (pi) for these initial densities. 
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Figure 9: Initial log-gamma densities Po(a, v, 1) from (18) are shown in the left panel 
for the uniform density r = v = 1, and also for r = 1 with v = 0.4 and 10; Figure [75| 
shows the evolution of the case with v — 10. The right panel shows the corresponding 



fractional decline with time of the population N(t)/N(fi) from (20) for these initial 



densities. Figure 10 shows the evolution from Po(a, 10, 1) 

— 100 




Figure 10: The evolution of probability density Pt(a, 10, 1) from (18) 
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Figure 11: Initial log-gamma densities Po(a, 7) are shown in the left panel for the 
uniform density t = V = 1, and also for r = 7 wif/i 1/ = 3, 10, 50 in graphs with 
increasing thickness. The right panel shows the corresponding fractional decline with 



time of the population N(t)/N(0) from (24) for these initial densities 
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3 Discussion 

These new results show how evolution of an inhomogeneous rate process reacts to 
perturbations of the uniform density using the smooth two parameter (is, r) family 
of log-gamma probability density functions [2] to represent unfitness measured by 
parameter a £ [0, 1]. We have illustrated the evolution from three initial states having 
central mean Eq(cl) — \ but different variances (cf. Figures [2] and [Zj) : 

initial uniformity (r = 1 with variance o"g(a) = A « 0.083) 

initial higher variance (r = 0.289 with Oq(o) — 0.164) 

initial lower variance (r = 2.24 with 0-9 (a) = 0.045). 

The graphs of these three starting log-gamma densities in the left hand sides of Fig- 



6|7|8 show widely differing shapes; the higher variance arises from peaks at the 
extremities of unfitness and the lower variance case has a symmetric bell-like curve 
about the mean. The evolution of features from the initial uniform density from 
Karev [5] is — cf. equation ^ above, 

te~ at , . 1 1 N(t) 1-e"* 

P (a) = l^P t (a) = —, E t{ a) = - + —, = _. (23 ) 



To this we add our analytic results (18 1, (19) for the initial log-gamma cases with 



t = 1, which agrees with (231 at v = 1 and the generalized hypergeometric solutions 
(21), (22) along with corresponding versions for integer r > 3. Our series approxi- 
mation agrees with the analytic results for the early development. Figure [6] shows 
that in each of our cases the evolution of mean unfitness E t (a) does indeed follow the 
expected equation ([5]), which expresses Fisher's law of natural selection [JJ. Also, in 
Figure §the very slow early evolution of its variance of (a) reflects the corresponding 
lack of curvature in the mean. In fact, though the decline of the mean fitness does 
change with the variation in starting conditions, the net population evolution N(t) in 
Figure [7] seems surprisingly stable under quite considerable alterations in the shape 
and variance of the initial density of inhomogeneity in fitness. The rate process evi- 
dently has a strong smoothing effect when we begin from a central mean, even though 
the variances differ widely. 

We can however obtain considerable changes in the evolution if we depart from an 
initial central mean. We illustrate this using our analytic solutions to allow evolution 
over longer periods. For example, the initial log-gamma density for r = 1, equations 
(18), (19), (20), gives the examples shown in Figure[9j Here we see the effect of v which 



acts as a location parameter for the initial density. Figure [TCH shows the evolution of 



the initial case with v = 10 into the probability densities Pt(a, 10, 1) from (18). 

Also, the initial unimodular log-gamma probability densities for r = 7, Pt(a, v, 7), 
give this generalized hypergeometric solution 

= L (t) = 7-F7O, v, v, u\ H-l, H-l, v+1, u+1, u+1, u+1, v+1; -t) (24) 

and again v serves as a location parameter. Figure[TT]illustrates cases for v — 3, 10, 50 
as well as that for an initial uniform density. 

Our approach has considered the determinate rate process for the evolution of 
individual types with an inhomogeneous fitness distribution; Baake and Georgii [3] 
consider the multi-dimensional systems of many types. Their mutation and differential 
reproduction processes yield n-dimensional systems of differential equations which 
they analyse using variational methods. 
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